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Abstract 

Background: Sex allocation of offspring in mammals is usually considered as a matter of chance, being dependent 
on whether an X- or a Y-chromosome-bearing spermatozoon reaches the oocyte first. Here we investigated the 
alternative possibility, namely that the oviducts can recognise X- and Y- spermatozoa, and may thus be able to bias 
the offspring sex ratio. 

Results: By introducing X- or Y-sperm populations into the two separate oviducts of single female pigs using bilateral 
laparoscopic insemination we found that the spermatozoa did indeed elicit sex-specific transcriptomic responses. 
IVlicroarray analysis revealed that 501 were consistently altered (P-value < 0.05) in the oviduct in the presence of 
Y-chromosome-bearing spermatozoa compared to the presence of X-chromosome-bearing spermatozoa. From these 
501 transcripts, 271 transcripts (54.1%) were down-regulated and 230 transcripts (45.9%) were up-regulated when the 
Y- chromosome-bearing spermatozoa was present in the oviduct. Our data showed that local immune responses 
specific to each sperm type were elicited within the oviduct. In addition, either type of spermatozoa elicits 
sex-specific signal transduction signalling by oviductal cells. 

Conclusions: Our data suggest that the oviduct functions as a biological sensor that screens the spermatozoon, 
and then responds by modifying the oviductal environment. We hypothesize that there might exist a gender 
biasing mechanism controlled by the female. 
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Background 

For many years gender allocation of offspring in mammals, 
including humans, has been regarded as a matter of 
chance, depending on whether an X- or a Y- chromosome- 
bearing spermatozoon reaches the oocyte first. Since an 
equal number of X- and Y- spermatozoa are produced dur- 
ing spermatogenesis [1], and fertilization is a random 
event, it stands to reason that in each generation equal 
numbers of males and females should be born. Evidence 
from the field and laboratory challenges this classic dogma 
and suggests that some kind of adaptive control of 
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offspring gender may exist in mammals [2,3] . Evidence for 
this ability exists in many invertebrates and some avian 
species are able to adjust their progeny sex ratio pre- 
dictably in response to environmental conditions [4]. 
Numerous factors such as population density, resource 
availability (famine), season, mother's age, levels of 
hormones, time of insemination and stress are known 
to influence the sex ratio in mammals [5-8]. However, 
the biological mechanism(s) through which mammals 
can bias the offspring ratio is still unknown. 

Several hypothetical mechanisms have been proposed 
to explain sex ratio skewing in mammals. On the male 
side, any shift from the expected 1:1 sex ratio among off- 
spring has been related to intrinsic differences in sperm 
motility, viability and fertilization ability of the two 
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gametes types [9]. On the female side, the condition of 
the reproductive tract and the penetrabiUty of the oo- 
cyte's zona pellucida, which varies according to the tim- 
ing of insemination relative to ovulation, have been 
suggested to influence differentially the ability of X- or Y- 
spermatozoa to fertilize oocytes [10]. Once fertilization 
has occurred, the milieu of the oviduct, and subsequently 
the environment of the uterus, may affect the develop- 
mental rates of XX-embryos and XY-embryos [11,12]. 

Given that the female investment in the offspring is 
considerably larger than that of the male, it is more 
probable that a mechanism to bias the offspring sex ratio 
is operated by the mother [13]. Furthermore, it is more 
likely that such a hypothetical mechanism would operate 
in the oviduct immediately before, or at the time of 
fertilization because it is less costly to females than other 
suggested mechanisms acting later during pregnancy 
[2,11]. The mammalian oviduct is the venue for import- 
ant reproductive events such as sperm and oocyte 
transport, sperm binding and release, fertilization and 
early-embryonic development [14]. In addition, the ovi- 
duct is implicated in the selection of spermatozoa, being 
capable of distinguishing between good and poor sperm 
quality [15]. 

Here, we address the question of whether the female may 
distinguish between the presence of X and Y-spermatozoa 
in the oviduct before fertilization occurs. We tested this 
possibility by examining whether the presence of the X- 
and Y-spermatozoa elicit different transcriptomic responses 
within the oviduct. To test our hypothesis we used an 



in vivo pig model that directly compared the oviduct con- 
taining Y-spermatozoa to the contralateral oviduct from 
the same animal, but containing X-spermatozoa (Figure 1). 
The advantages of this model were: 1) that minimize the 
confounding factors known to bias the sex ratio [2] since 
both oviducts analyzed were from the same animal and 
therefore were under the same nutritional, health and hor- 
monal environment, and 2) that avoid the possibility that 
oocytes could mask the oviductal responses towards X- 
and Y-spermatozoa, because only sows showing multiple 
pre-ovulatory follicles were selected for this study. It has 
previously been reported that, like spermatozoa, oocytes 
elicit distinct proteomic alterations [16]. 

Our study add a complete new layer of competition to 
the mating game, since up to date most studies of off- 
spring sex ratio are based on epidemiological studies, 
showing a traditional maternal dominance or lately a 
male influence in specific species. We open up a new 
perspective in "the battle of the sexes", suggesting that 
this battle starts in the oviduct and providing the first 
molecular evidence of a sex-specific sperm recognition 
system in the oviduct. 



Results and discussion 

Our work showed that the presence of X- and Y- 
spermatozoa did indeed elicit different transcriptomic 
responses within the oviduct (Figure 2A). Around 2% of 
transcripts (501 out of 24123 probes from Affymetrix 
Porcine Chip) were consistently altered (P-value < 0.05) 
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Figure 1 Schematic representation of the experimental design. Sows were subjected to laparoscopic surgery. To prevent X- and 
Y-spermatozoa migration between oviducts, both uterine horns were cut using titanium staples. Then, one oviduct was inseminated with 
X-spermatozoa and the contralateral oviduct was inseminated with Y-spermatozoa (3x10^ spermatozoa/100 (il) from the same animal. Twenty-four 
hours following laparoscopic insemination, oviductal tissues containing X- and Y-sperm samples were collected from each side of the 
reproductive tract in all animals. 
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Figure 2 The presence of Y-spermatozoa elicited different transcriptome response within the oviduct when compared to X-spermatozoa. 

A: Cluster heat map analysis of the transcriptional profiles obtained from oviductal samples inseminated with X-spermatozoa and Y-spermatozoa. 
Each row represents a different gene, and each column displays gene expressions at different samples (XI -X4 for oviductal samples inseminated 
with X-spermatozoa; Y1-Y2 for oviductal samples inseminated with Y-spermatozoa). Data values displayed as yellow and blue represent elevated 
and reduced expression, respectively. B: A Volcano plot depicting significant changes in gene expression between oviductal samples inseminated 
with Y-spermatozoa and X-spermatozoa. Each of the 23,124-oligonucleotide probes is represented by a dot in the graph. The x-axis represents the fold 
change and the y-axis represents the statistical significance (-loglO of p-value). Transcripts showing significant differences in gene expression 
(501 probes, p < 0.05) are above the broken line. 



in the oviduct in the presence of Y-chromosome-bearing 
spermatozoa compared to the presence of X-chromosome- 
bearing spermatozoa (Figure 2B). From these 501 transcripts, 
271 transcripts (54.1%) were down-regulated and 230 
transcripts (45.9%) were up-regulated when the Y- 
chromosome-bearing spermatozoa was present in the 
oviduct. A complete list of the transcripts altered in the 
oviduct inseminated with Y- chromosome bearing sperm- 
atozoa compared to X -chromosome bearing spermatozoa 
is presented in the Additional file 1. 

To obtain a biologically meaningful overview of the al- 
tered transcripts in the presence of Y-chromosome-bearing 
spermatozoa compared to X-chromosome-bearing sperm- 
atozoa, genes differently expressed were organized into 
different categories and subcategories according to KEGG 
database hierarchy. The functional categories with higher 
number of genes were: signal transduction, immune sys- 
tem, digestive system and endocrine system. The pathways 
in which these altered transcripts were involved with 
are presented in Table 1. Interestingly, a higher number 
of genes involved in signal transduction and immune 
system were up-regulated (60-70%) in the presence of 
Y-chromosome bearing spermatozoa when compared to 



X- chromosome-bearing spermatozoa. Other interest- 
ing subcategories with high numbers of transcripts in- 
volved were: nervous system, cell growth and death, cell 
communication, signalling molecules and interaction, 
folding, sorting and degradation and transport and ca- 
tabolism. The results of all data pathways classification 
are available in Additional file 2. 

Our data provide the first evidence to show how sperm- 
atozoa carrying the Y- or X-chromosome can modulate 
the oviductal response by activating specific signalling 
pathways in a gender specific manner. These data imply 
that the female reproductive tract recognizes the presence 
of X- or Y-chromosome-bearing spermatozoa in the 
oviduct before fertilization occurs. On this basis, we 
hypothesize that a sex-specific sperm recognition sys- 
tem exists in the female reproductive tract. Further- 
more, we propose that this sperm recognition system 
can modulate the gender selection of the offspring. Our 
hypothesis challenges two long-held assumptions: 1) on 
the female side, that the reproductive tract is a passive 
participant in sperm selection [17] and 2) on the male 
side, that there are no differences in morphology, meta- 
bolic activity or functional abilities of X and Y sperm that 
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Table 1 Pathways and altered genes Involved in signal transduction, immune system, digestive system and 
endocrine system 



KEGG category KEGG KEGG pathways Transcripts 

subcategory 



Environmental Information Signal 


Wnt signaling pathway 


PPP2R5C, FBXW1 1, SMAD4, WNT5B, PRKCA 




MARK signaling patiiway 


TGFB1, FGF12, CACNA1B, CRK, FGFR1, PRKCA, MAP4K3, 
MAP3K1 1 




TGF-beta signaling pathway 


TGFBl, SMAD5, LTBP1, SMAD4, BIV1P4 




mTOR signaling pathway 


ULK2, PRKAA1, HIF1A 




Calcium signaling pathway 


PDE1A, CACNAIB, GNAL, HTR2A, ERBB4, ADCY9, PHKAl, 
PnPlR PRKCA l\inS7A 




FrhR t;Inna inn piaTn\A/a\/ 
Liuu oiyi lain ly ua li i vvay 


CRK FRRR4 CDKNIR PRKCA 




Pnnt;)~inatirl\/ inr\c:itin t;inna \r\n c\/c:tprn 
r 1 luiui laLiuyiM iujilui jiyi icim ly jyoLCi i i 


nCKI PRKCA INPP4R 




jai\ JirAi jiyi iciiii ly ud LI 1 vvci y 


LEPR, LIF 




1— IpHnphriD <;inn3!inri ri3l'h\A/;^\/ 
1 icuyci ivJy jiyi iciiii ly uoli ivvay 


FRXW1 1 WNTSR RMP4 




V cvjr jiy I idi 1 1 ly [Jd LI 1 vvdy 


PRKCA 


<yi ijui iijiiiui -jyj ici i ij 1 1 1 1 1 1 1 ui ic 

System 


1 1 1 LCroLi iidi iiiiiiiLJiic iicrLVVUii\ ilji ly/* 

production 


TGFBl 




Chemokine signaling pathway 


CCL8, CRK, ADCY9, VAV2 




Hematopoietic cell lineage 


CD55 




Complement and coagulation cascades 


CD55, TFPI, C5AR1 




Fc gamma R-mediated phagocytosis 


CRK, PPAP2A, PRKCA, VAV2 




Leukocyte transendothelial migration 


CTNNA3, CLDNIO, IV1MP9, PRKCA, VAV2 




Natural killer cell mediated cytotoxicity 


FCERIG, PRKCA, VAV2 




Fc epsilon Rl signaling pathway 


FCERIG, BTK PRKCA, VAV2 




Toll-like receptor signaling pathway 


MYD88, IRF7 




B cell receptor signaling pathway 


BTK, VAV2 




RIG-l-like receptor signaling pathway 


ATG5, IRF7 




Cytosolic DNA-sensing pathway 


IRF7 




T cell receptor signaling pathway 


VAV2 


Digestive 


Vitamin digestion and absorption 


SLC23A1 


System 


Mineral absorption 


CYBRD1, ATP1B2, ATP1A2, ATP7A 




Fat digestion and absorption 


PPAP2A 




Salivary secretion 


ATP1B2, ATP1A2, ADCY9, PRKCA 




Gastric acid secretion 


ATP1B2, ATP1A2, SSTR2, ADCY9, PRKCA 




Pancreatic secretion 


ATP1B2, ATP1A2, ADCY9, CLCAl, PRKCA 




Bile secretion 


ATP1 B2, ATPl A2, ABCBl 1, ADCY9, ABCBl 




Carbohydrate digestion and absorption 


ATP1B2, ATP1A2, SLC2A2 




Protein digestion and absorption 


ATPl B2, ATPl A2, C0L11A1 


Endocrine 


Progesterone-mediated oocyte maturation 


CDC27, SPDYA, ADCY9, CCNAl 


System 


Insulin signaling pathway 


CRK, PCK2, PRKAAl, PHKAl, IRS4, FL0T2 




Adipocytokine signaling pathway 


ACSBGl, PCK2, LEPR ADIPOQ, PRKAAl, IRS4 




PPAR signaling pathway 


ACSBGl, PCK2, ADIPOQ, FABP7, PLTP 




GnRH signaling pathway 


ADCY9, PRKCA 




Melanogenesis 


ADCY9, WNT5B, PRKCA 



could be used by the female tract to recognize the two 
types of sperm [18]. 

Recent research has tended to suggest that the female 
tract is an active participant in the sperm selection 
process [15,19]. Furthermore, the oviduct is emerging as 



a chief protagonist in the sperm selection process, facili- 
tating/or inhibiting sperm transport and allowing only se- 
lected sperm to reach the oocyte. This poses the question 
of whether spermatozoa are pilots or merely passengers in 
the sperm journey to the oocyte. It is estimated that from 
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the approximately 30 billions of spermatozoa deposited 
into the female reproductive tract after insemination, only 
around 1000-5000 spermatozoa enter the oviduct in the 
pig and bind to the oviductal cell surface [20] . Once in the 
oviduct the spermatozoa are subjected to several further 
selection processes before being able to interact with oo- 
cytes. This oviductal selection process seems to be based 
on the intrinsic integrity and information content of the 
sperm DNA or/and based on more subtle properties that 
reflect the individual spermatozoon [15]. Thus, we specu- 
late that the oviduct could also differentiate between X 
and Y-spermatozoa. 

Researchers have for many years investigated the differ- 
ences in size, shape, motility or differentially expressed pro- 
teins between X- and Y-spermatozoa [21-23] and no 
conclusive findings have emerged, aside from the very small 
2.8-4% difference in genetic material (depending on spe- 
cies). This small difference in DNA content, due to the arm 
of the X chromosome that is not present on the Y, is tightiy 
packaged in a semi-crystalline form inside the spermato- 
zoon [24] and does not seem to be accessible for evaluation 
by external systems in the laboratory. Our new findings 
open the possibility that the oviductal recognition of X- and 
Y-spermatozoa may be intimately related to currently un- 
known differences in morphology or metabolism of X- ver- 
sus Y-bearing sperm. Recently researchers have identified 
different topographic characteristics on the head of X- and 
Y-spermatozoa, viewed on a nanometric scale using atomic 
force microscopy [25]. Another study has pointed to differ- 
entially expressed proteins between bull X- and Y- 
spermatozoa involved in energy metabolism, stress resist- 
ance and cell defence [23]. Chen and colleagues [23] even 
suggested that there could be differences in the way that 
energy is produced and varied vulnerabilities to environ- 
mental changes between X- and Y-spermatozoa. It is 
known, for example that different metabolic rate is related 
to dissimilar production of ROS [26]. Thus, we speculate 
that the different sperm metabolism identified by Chen and 
colleagues [23] could lead to the release of distinct amount 
of ROS substances or diverse ROS metabolites to the ovi- 
duct by each type of spermatozoa. To investigate this we 
measured sperm motOity between X and Y-sperm samples 
(89.94 ± 1.26% and 91.55 ± 0.78%, respectively), sperm viabil- 
ity (89.57 ± 1.3 and 90.57 ± 1.61 respectively) and intracellu- 
lar ROS generation (124.74 ± 17.96 and 126.89 ± 19.54 FU/ 
lO'^^ live spermatozoa, respectively). As our analyses did not 
reveal any differences we conclude that the oviducts ability 
to differentiate between the two types of spermatozoa must 
be based on other mechanisms. Nevertheless, our data dem- 
onstrated that the specific oviductal response to each type of 
sperm was not as a result of a higher number of dead sperm 
in one type of sperm sample compared to the other or a 
higher percentage of sperm damage by the flow sorting 
method and therefore a higher production of ROS. 



Immunological strategies to gender selection have also 
been proposed since the finding of a family of gene 
products encoded in, or controlled by, the Y chromo- 
some that is only present on male's cell surfaces [27]. 
This idea, that sex selection in mammals occurs through 
a specific immune response, is in line with our results. 
Bioinformatics analysis of our microarray data showed 
that genes involved in signal transduction and immune- 
related genes were the most representative of the altered 
genes in this study (Figure 3). Moreover, a higher number 
of genes involved in these pathways were up-regulated 
(60-70%) in the presence of Y spermatozoa when com- 
pared to X spermatozoa (Supplementary data). Particu- 
larly interesting was the fact that the immune system 
seems to be involved in the recognition of X and Y sperm- 
atozoa within the oviduct. It is logical to assume that the 
immune system plays a key role in ensuring tolerance to 
spermatozoa in the maternal tract. Under normal circum- 
stances, when the maternal tract is exposed to pathogens 
or a non-self-entity, the immune system responds in an 
aggressive manner. However the spermatozoon, which is a 
non-pathogenic "foreign invader" to be destroyed (non- 
self-entity), is accepted in the maternal tract and is guided 
to the oocyte. We suggest that the immune system might 
act as a molecular screening process in the oviduct that al- 
lows only preferred X- or Y-spermatozoa to reach the oo- 
cyte. Interferon regulatory factor 7 (IRF7) and Chemokine 
(C-C motif) ligand 8 (CCL8), which are involved in the 
immune system as signalling molecules were further con- 
firmed using qPCR (Figure 4). IRF7 is associated with spe- 
cific families of pattern recognition receptors such as Toll- 
lilce receptors (TLRs), RIG-I-lilce receptors (RLRs) and 
Cytosolic DNA-sensing pathway [28]. All of these path- 
ways are responsible for detecting microbial pathogens or 
foreign DNA from invading microbes or host cells and 
may also participate in the recognition between X and Y 
sperm. The current data further corroborated and con- 
firmed our previous report showing a distinct response of 
oviduct to spermatozoa of unsorted ejaculates [16]. How- 
ever, due to the limited sample size, the results presented 
here should be interpreted with caution. Our future exper- 
iments would be directed towards increasing the sample 
size of the study as well as establishing an in vitro based 
system for understanding the mechanism(s) mediating 
this process. 

Conclusions 

The present investigation demonstrates for the first time, 
distinct alterations of oviductal gene expression in response 
to X and Y spermatozoa. These changes imply that the ovi- 
duct is able to distinguish between X and Y-spermatozoa 
and thereby fine-tune its physiology and gene expression 
in response. We propose that the oviduct functions as a 
biological sensor that screens the spermatozoa, allowing 
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Figure 3 Transcripts modulated by X-spermatozoa and Y-spermatozoa in the oviduct organized Into functional categories. Transcripts 
differentially expressed in oviductal samples inseminated with X-spermatozoa compared to Y-spermatozoa organized into functional categories 
on the basis KEGG PATHWAY database. 



only a preferred cohort to proceed. Our findings point to 
the existence of a sex-specific-sperm recognition system in 
the oviduct that alerts the mother to the presence of the X 
or Y-spermatozoa, thereby expanding currently proposed 
hypotheses that influence sex ratio populations. 

Methods 

Experimental design 

To determine the alterations in the oviduct transcriptome 
in response to X and Y sperm our unique experimental 
model together with Affymetrix microarray technology 
were used. Figure 1 provides a schematic overview of the 
experimental model employed. Four biological replicates 
were performed (n = 4 sows) and a total of 8 arrays were 
used for microarray study (4 arrays for oviduct samples 
containing Y- bearing spermatozoa and 4 arrays for ovi- 
duct samples containing X-bearing spermatozoa). To 
validate the microarrays results, two transcripts were 



selected (Chemokine (C-C motif) ligand 8 (CCL8) and 
Interferon regulatory factor 7 (IRF7)) and gene expres- 
sion was analyzed by qPCR. For this experiment, 4 bio- 
logical replicates (n = 4 sows) and 3 technical replicates 
were performed. 

Animals 

Weaned crossbred sows from Landrace x Large White 
genetic line (from two to six parities) were selected for 
these experiments. All experiments were performed after 
obtaining approval from the Ethical Committee for 
Experimentation with Animals of the University of Murcia, 
Spain (385/2008). 

Detection of oestrus and ovarian status 

Oestrus detection was carried out once a day, 2 days 
after weaning, by exposing females to a mature boar and 
applying manual back-pressure. Females that showed a 
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Figure 4 Validation of the microarray results by qPCR analysis. 

CCL8 (chemokine (C-C motif) ligand 8) and IRF7 (interferon 
regulatory factor 7) expression values (normalized based on l]-actin 
and Ubiquitin B expression values) in oviductal samples inseminated 
with X-spermatozoa compared to oviductal samples inseminated 
with Y-spermatozoa. The expression of both transcripts in the 
oviductal samples inseminated with Y-spermatozoa was significantly 
different from that of the oviductal samples inseminated with 
Y-spermatozoa (P < 0.05). 



Standing estrous reflex were considered to be in heat 
and the ovaries were scanned by transrectal ultrasonog- 
raphy [29]. Only sows showing multiple pre-ovulatory 
follicles (diameter of antrum >6 mm) were selected for 
experiments. Inseminations were carried out within 2-3 h 
after the ultrasonography and oviduct samples collection 
took place 24 h later (approximately 12-24 h before ovula- 
tion took place in most of the sows). 

Semen collection and flow cytometric sorting 

Sperm-rich fractions from three fertile mature boars, 
that had previously sired offspring, were collected by 
gloved-hand method and extended in Beltsville Thawing 
Solution (BTS) to 150 x 10*^ spermatozoa/mL [30]. After 
collection, samples were evaluated for normality (motil- 
ity >80%, membrane integrity > 85%, total sperm count 
per ejaculate > 20 x 10^; acrosomal abnormalities < 10%, 
abnormal sperm morphology < 15%) [31]. 

The extended semen was then processed for sperm 
sorting following the general procedure as previously de- 
scribed [32] . Briefly, 1 mL of extended semen was stained 



with Hoechst-33342 fluorophore (0.3 mM per 1 x 106 
spermatozoa) and incubated for 1 h in darkness at 35°C. 
After incubation, samples were passed through a 30 mm 
nylon mesh filter to remove debris or clumped sperm- 
atozoa. Then, 1 [^l/ml of food dye (1 mg/ml; Warner 
Jenkinson Company Inc., St. Louis, MO, USA FD&C 40 
solution) was added to the sample to identify dead 
sperm and kept 15 minutes at room temperature in 
dark. The stained spermatozoa were separated using a 
high speed MoFlo SX flow cytometer/sperm sorter 
(Dako Colorado Inc., Fort Collins, Co, USA), equipped 
with a solid state laser (ultraviolet wavelength, 351- 
363 nm) at 175 mW (Spectra Physics 1330, Terra Bella 
Avenue, Mountain View, CA, USA). Spermatozoa were 
separated into X- and Y- chromosome bearing popula- 
tions with 91% purity of Y-spermatozoa and 93% of X- 
spermatozoa. Both the X sperm and the Y sperm popula- 
tions were recovered from the three boars. Consequently, 
in one side the three X-insemination doses were pooled 
and in the other side the three Y-insemination doses from 
the three boars were also pooled and both doses finally di- 
luted in BTS to 3 X 10^ spermatozoa/ 100 |il before laparo- 
scopic insemination. 



Sperm motility, viability and functionality of X- and 
Y-sperm samples 

The motility of X- and Y-sperm samples was evaluated ob- 
jectively using a computer-assisted analysis system (ISAS; 
Proiser R -i- D, Paterna, Spain). The sperm motility vari- 
able recorded was the overall percentage of motile sperm- 
atozoa (average path velocity (VAP) = 20 mm/s) [31]. 

The viability of the spermatozoa was evaluated by sim- 
ultaneous cytometric assessment of the plasma and acro- 
somal membrane integrity by using a triple-fluorescence 
procedure as described by Martinez- Alborcia et al., [33]. 
Briefly, a 100 |^1 sperm sample (30 x 106 cells/ml in 
PBS) was transferred to culture tubes containing 2.5 [A 
H-42 (0.05 mg/ml in PBS), 2 ^il propidium iodide (PI, 
0.5 mg/ml in PBS, Molecular Probes Europe BV, 
Leiden, The Netherlands) and 2 i^l fluorescein-conjugated 
peanut agglutinin (PNA-FITC, 200 \xg/m\ in PBS). The 
samples were mixed and incubated at 38°C in the dark for 
10 min. Immediately before analysis by flow cytometry, 
400 |il PBS was added to each X- and Y- sperm samples 
and mixed. The fluorescence spectra of PI and PNA-FITC 
was detected using a 670 nm long-pass (LP) filter and a 
530/30 nm BP filter, respectively. The spermatozoa ana- 
lyzed were categorized into four categories: (1) intact 
plasma and acrosomal membranes (PI-/PNA-FITC-); 
(2) intact plasma membrane and damaged acrosome 
(PI-/PNA-FITC-H); (3) damaged plasma membrane and 
intact acrosome (PI-I-/PNA-FITC-); or (4) damaged plasma 
and acrosomal membranes (PI-I-/PNA+). The viable 
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spermatozoa exhibited intact plasma and acrosomal mem- 
branes and this was expressed as a percentage of the total 
cells. 

The functionality of the spermatozoa was assessed by the 
intracellular ROS production in X and Y-sperm samples. 
The intracellular production of ROS by each sperm sample 
was measured using 5-(and- 6) chloromethyl-20,70-dichlor- 
odihydrofluorescein diacetate acetylester (CM-H2DCFDA), 
described by Martinez-Alborcia et al., [33]. 

CM-H2DCFDA is freely permeable across cell mem- 
branes and becomes incorporated into the hydrophobic 
regions of the cell. Upon entering the cell, the acetate 
moiety of CM-H2DCFDA is cleaved by cellular esterases 
to leave the impermeant and non-fluorescent molecule 
20,70-dichlorodihydro- fluorescein (H2DCF). The H2DCF 
is oxidized by hydrogen peroxide (H202) into dichloro- 
fluorescein (DCF), which fluoresces at 530 nm following 
excitation at 488 nm. For each X-and Y-sperm sample, 
two different 50 \il aliquots of mTBM-diluted spermato- 
zoa were diluted in 950 ^il PBS containing (1) 1.25 |il H-42 
(0.05 mg/ml in PBS), 1 ^il PI (1 mg/ml in PBS), 1 ^1 
H2DCFDA (1 mM in DMSO) and 1 ^1 tert-butylhydrogen 
peroxide (1 mM in purified water) to induce oxidative 
stress (first aliquot; used to measure induced ROS for- 
mation) or (2) 1.25 ^1 H-42, 1 ^il PI and 1 ^1 H2DCFDA 
(second aliquot; used to measure basal ROS formation). 
The samples were incubated at 38°C in the dark for 
10 min before flow cytometry. The mean fluorescence 
intensity of DCF (induced minus basal) was expressed 
as fluorescence units (FU) per 1012 live spermatozoa. 

Intraoviductal laparoscopic insemination 

Intraoviductal laparoscopic inseminations were carried 
out within 2-3 h of the ultrasonography. Sows were se- 
dated by azaperone administration (2 mg/kg body weight, 
i.m.), anaesthetized with sodium thiopenthal (7 mg/kg 
body weight, i.v.) and maintained under anesthesia with 
isofluorane (3.5-5%). Intraoviductal laparoscopic insemi- 
nations were carried as described by Alminana et al., [34] . 
Briefly, each sow was placed in the supine position and a 
pneumoperitoneum was established. The abdominal cavity 
was insufflated with C02 to 14 mmHg. Two accessory 
ports were placed in the right and left part of the hemi ab- 
domen, which provided access for laparoscopic Duval for- 
ceps for manipulating the uterine horn and grasping the 
oviduct for the insemination, respectively. To prevent X 
and Y sperm migration between oviducts, both uterine 
horns were cut using titanium staples (EndoGIA Universal 
60/4.8; Tyco heathcare, Mansfield, MA). To perform the 
insemination, the oviduct was grasped with the Duval for- 
ceps in the isthmus region and the sperm dose containing 
Y sorted spermatozoa (3 x 105 spermatozoa/100 i^l.) was 
flushed into one oviduct (above of the ampullar region in 
direction to isthmus). The procedure was then repeated in 



the contralateral oviduct but injecting the X sperm dose 
containing the X-chromosome bearing sperm (Figure 1). 
After both oviducts were inseminated, the trocars were re- 
moved and minor suture was required. Following surgery, 
sows were returned to their accustomed environment. 

Sample collection and RNA preparation 

Twenty-four hours following laparoscopic insemination, 
oviduct tissues containing X and Y sperm populations 
were collected from each side of the reproductive tract 
in each animal. Sows were sedated as previously described. 
To avoid the presence of oocytes that could mask the ovi- 
duct response to X and Y sperm, only oviduct samples 
from sows showing no signs of ovulation were collected. 
In addition, oviducts were flushed with Phosphate Buff- 
ered Saline medium (PBS, 30 ml) and the absence of oo- 
cytes was verified by careful examination of oviduct 
flushings under a stereomicroscope. In cases where it was 
unsure whether the sow had ovulated oviduct samples 
were discarded. 

Oviducts were opened longitudinally and epithelial cells 
were isolated by scraping the mucosal epithelial layer with 
a glass slide. Scraped cells from the uterine horn samples 
were transferred immediately to Tri Reagent (Sigma, 
Sigma-Aldrich Co, Madrid, Spain), homogenised, snap- 
frozen in liquid nitrogen and stored at -80°C until fur- 
ther processing. 

Total RNA was isolated using a standard procedure in- 
volving phenobchloroform extraction followed by ethanol 
precipitation. The quantity (NanoDrop 1000 spectro- 
photometer) and the quality (Agilent 2100 Bioanalyser; 
Agilent Technologies) of the RNA samples were ana- 
lysed. Only samples with satisfactory quality as indi- 
cated by the absence of degradation of the ribosomal 
RNA were used for microarrays and quantitative Real- 
Time Reverse Transcriptase-Polymerase Chain Reaction 
(qPCR). 

Microarrays hybridization 

Affymetrix Porcine Genome gene expression arrays 
(Affymetrix, Santa Clara, CA) were used in this study. 
Total RNA samples were prepared according to the Affy- 
metrix Technical Manual (www.affymetrtx.com). Briefly, 
200 ng of total RNA was converted into cDNA using an 
oligo(dT) which also carries the binding site for T7 RNA 
polymerase. Superscript II (Affymetrix) was used to carry 
out this reaction. After first strand synthesis, residual 
RNA was degraded by addition of RNaseH and a double- 
stranded cDNA molecule was generated using DNA Poly- 
merase I and DNA ligase. This double stranded molecule 
was used as a substrate for the T7 RNA polymerase to 
produce multiple copies of the cRNA using the Affymetrix 
IVT labelling system. The cRNA molecules produced in- 
corporated biotin labelled ribonucleotides, which acted as 
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a target for the subsequent detection of hybridization, 
using fluorescently labelled streptavidin. 13 i^g of cRNA 
molecules were heat fragmented and injected to the Por- 
cine GeneChips in a hybridization solution according to 
the Affymetrix protocol. Hybridization took place over- 
night in a rotating hybridization oven at 60 rpm, 45°C for 
16 hours. The GeneChip arrays were washed using the 
Affymetrix Fluidics Station. After washing, the GeneChip 
arrays were scanned using Affymetrix GC3000 scanner. 
The resultant images were analysed using the Microarray 
Suite software version 5.1 (Affymetrix). At the detection 
level each probe set was designated as present, absent or 
marginal. Only present transcripts were considered 
expressed. Microarray experiments were carried out ac- 
cording to MIAME guidelines and the complete experi- 
mental data can be obtained online from the NCBI Gene 
Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) 
submission number GSE47139. 

Microarray data and bioinformatics analysis 

Microarrays analysis was performed using Taverna work- 
flow management system (http://www.taverna.org.uk.) 
[35,36]. Microarrays data were normalized using Robust 
Multichip Average (RMA) method [37]. Differentially 
expressed genes between oviduct inseminated with X 
spermatozoa and Y spermatozoa were detected through 
t-tests that were applied to normally distributed data 
using the Limma R/Bioconductor packages [38] with 
Benjamini-Hochberg false discovery rate and multiple 
testing correction to control for Type I errors [39]. From 
the list of differentially expressed genes, the genes that 
pass the specified threshold for p-value = 0.05 were se- 
lected. Using Affymetrix Porcine Annotation in combin- 
ation with the improved annotation provided by Tsai 
et al., [40] for these porcine arrays, differentially expressed 
transcripts from our microarray data were annotated. To 
obtain a biologically meaningful overview of the altered 
transcripts in the presence of Y- chromosome-bearing 
spermatozoa compared to X-chromosome-bearing sperm- 
atozoa, genes differently expressed were organized into 



different categories and subcategories according to KEGG 
database hierarchy. 

Quantitative real-time reverse transcriptase polymerase 
chain reaction 

Gene expression profiles derived from microarray analyses 
were confirmed using quantitative real-time reverse tran- 
scriptase polymerase chain reaction (qPCR). The primers 
used for qPCR are listed in Table 2. Amplified PGR prod- 
ucts were sequenced with forward and reverse primers to 
verify the resulting product. 

Total RNA from the oviduct samples (inseminated 
with Y sperm or inseminated with X sperm) was treated 
three times with DNase I (DNA-free kit; Ambion.) to re- 
move genomic DNA contamination from samples. First- 
strand cDNA synthesis was performed using High Capacity 
cDNA Reverse Transcription Kit (Appied Biosystems). 
Negative controls were prepared without inclusion of the 
enzyme (non-reverse transcription controls, RT controls). 
Reverse transcriptase PGR (RT-PCR) was performed ac- 
cording manufacturer instructions. To evaluate the size of 
the PGR products 10 |il of each sample was resolved on a 
1.2% agarose gel and electrophoresis was performed with 
Ix TAE buffer and a voltage of 110 V for 40-50 min. The 
bands were visualized by using an ultraviolet transillumin- 
ation, and digital images were obtained. 

SYBR Green Jump Start (Sigma) master mix (containing 
10 |il SYBR Green, 7 \il H20, 1 |il of forward and reverse 
primers and 1 |il cDNA) was added to each well of PGR 
plate and amplification was performed under the follow- 
ing conditions: 40 cycles of 95° for 30 s, 55° for 1 min and 
72° for 1 min. Samples without template and RT controls 
(without the addition of enzyme) for each primer set were 
included to identify contamination. Triplicate measure- 
ments for each group of samples were carried out. 
Quantitative PGR was performed using Mx3005P QPGR 
(Stratagene, Waldbronn, Germany). The quantification 
data were analyzed using MxPro QPCR software version 
4.01. Quantitative PCR data were analyzed using the com- 
parative CT method [41]. The results were expressed as 
mean ± SEM arbitrary gene expression values, normalized 



Table 2 Primers used for qPCR analysis 

Gene symbol Affymetrix porcine probe Primer Sequence Product size (pb) 

CCL8 SSc.9957.1Al_at Forward 5' GCGAGATGGCAnTCTCTCT 3' 119 

Reverse 5' CACACTTCGGCTTACAAGAGG3' 

IRF7 S5c.2573.9.1_51_at Forward 5' GCTGGATGAAGCCAGAACA 3' 97 

Reverse 5' GGCCCAGGCCTTAAAGAT 3' 

Ubiquitin B Reference gene Forward 5' GTCTGAGGGGTGGCTGCTAA 3' 85 

Reverse 5' TGGGGCAAATGGCTAGAGTG 3' 

P-actin Reference gene Forward 5' CCTCCCrGGAGAAGAGCTA 3' 131 

Reverse 5' CTTCATGATGGAGTTGAAGGT 3' 
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on the basis of the two reference expression (fi-actin and 
Ubiquitin B). Statistical analysis was performed using 
paired T-test to evaluate the significance of difference be- 
tween expression values of oviduct inseminated with Y 
sperm versus oviduct inseminated with X sperm (in SPSS, 
version 14.0 (SPSS Inc., Chicago, IL)). The threshold for 
significance was set at p < 0.05. 

Availability of additional files 

The data set supporting the results of this article is avail- 
able in the additional files and also is available online 
from the NCBI Gene Expression Omnibus (http://www. 
ncbi.nlm.nih.gov/geo/) submission number GSE47139. 
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